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FACILITY: MATH LIBRARY 
++ 
ABSTRACT: 
MTHSEXP is a function which returns the single floating point 


exponential of its single pres tacos floating point argument. 
The call is standard call-by-reference. 


VERSION: 0 
HISTORY: 
AUTHOR: 
Peter Yuo, 15-Oct-76: Version 0 
MODIFIED BY: 


O-1 Peter Yuo, 22-May-77 
1-012 Jeffrey C. Wiener 22-Feb-82 
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VXSEXP ; Single Precision Floating Ex one BS SEP-1984 4: AX/VMS Macro v04-00 P 
wt fin] ORY ; Detailed Current Edit Hi History <* -SEP-1984 1 9% 3 3 MTHRTL.SR eu UVXEXP.MAR; 1 — 3) 


1 -SBTTL HISTORY ; Detailed Current Edit History 
; Edit History for Version 0 of MTHSEXPDEXP | 
5 ; 0-1 Code saving after code review mer h 1977 
§ ; O-2 Finish error handling 10-June-1 
a0 ; 3 change RET to rss in ERROR; fix undefined GLOBLs | 
0 ; 0-6 MINSSERROR changed to to M+ saeetchatateat 
00 60 ; MTH$_... change | 
00 61 ; Changed error oth fy mechanism. Put error result in RO before 
Boh 66 3 calling MTHS$SIGNAL in order to allow user mggi ty ercor result. 
0 65 ; o-2 Declared PSECTs and use SF$W_SAVE_PSW. TNH 20-Dec-77 
009 be ; 0-8 Invoke SSFDEF. TNH 20-Dec-77 
0990 6 ; Edit History for Version 1 of MTHSEXP 
0000 68 ; 1-1 Split single and double precision routines into two parts; 
0000 +4 ; Used gore ogcurate and faster algorithms provided by M. Payne. 
0000 0; 7 23-Jan- 
0000 71 ; 1-3 Fixed bug causing unexpected integer overflow. JMT_24-MAR-78 
0000 if 3 1-46 Fixes for better accurac Use ADDD. JMT 11-Apr- 
0000 73; 1-5 Move .ENTRY mask definition to module header. TNH 14-Aug-78 
4444 74 ; 1-006 - peeete version number and co _* ht notice. JBS _— -78 
000 75 ; 1-007 - rie symbols MTH__FLOUNDMA MTH__FLOOVEMAT t 
44 4 3 “at FLOUNDRAT and MTHS$K "FLOOVEMAT, “Sospoebieniine 
$000 78 ; 1-008 - aad 5.0% ¥ “abe PSECT directive. JBS 22-DEC-78 | 
0000 79 ; 1-009 = Declare externals. SBL 17-May-1979 
0000 80 ; 1-010 = Included logic to avoid the loss in si nificance in EMOD for 
44 4 3 arguments greater than 2**4, RNH 23- 
444 ¢ 3 1-011 - chang ~ ‘dl 12, s* in calls to MTHS$SIGNAL and MTH$S$JACKET_TST 
; -Se 
0000 54 : 1-012 - changed 4s jaterehiae to D_floating point operations to 
0000 85 ; floating point operations. This was done to conform to 
0000 86 ; the policy that G_floating, not D_floating, should be used 
0000 rH 3 to back-up F g/ boat ing pe point ggreretions. I also corrected 
0000 8 ; some comments. JCW 22-F | 
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; Declarative Part of Modu 


-DSABL 
-E 
E 
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; MACROS: 
SSFDEF 


-PSECT 


OWN STORAGE: 
CONSTANTS: 


TABLE: 
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INCLUDE FILES: 
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; Declarative Part of Module 


MTHJACKET.MAR 


EXTERNAL SYMBOLS: 


; PSECT DECLARATIONS: 


_MTHSCODE 


none 


; Table to be used for scalin 
; verified by M. Payne using 
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GBL 
XTRN MTHSK_FLOUNDMAT 
eEXTRN MTHSK-FLOOVEMAT 
EX MTHSSSIGNAL 
-EXTRN MTHSSJACKET_TST 


EQUATED SYMBOLS: 
ACHASK = “MCV, R2 R3, RQ> ss; 
Ln. 8 “051 


her 


; SIGNAL SEVERE error 


fegtster saving mask and IV enable 
; Extension for operand in EMODF 


: define SF$ (stack frame) symbols 
PIC,SHR,LONG, EXE ,NOWRT 


3; program section for math routines 


These constants here have been 
program ROOT16 on PDP-10. 


: 2*(0/16) = 1.0 
> : Senta/i6y 
3 ote 
‘ 3 na(8/18) 
> 2*2(4/16) 
: 2**(5/16) 
3 mets ADS 
GE SG 
s an 
& ; ne(9/16) 
2 > 2*2(10/16) 
4 3: 2e*(11/16) 
; 3 eb TAY 
3 2**(15/16) 
7 > 2e*(14/16) 
5 3; 2e*(15/16) 


This was used when D_float backed up F_floating 
Table to be used for scaling. These constants here have been 
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: Singl 4 02:04:59 VAX/VMS Macro v04-00 
beCLARAT! ; Declarative Part of Modu gS Ei obec 99:08 8e MTHRTL.SRCJUVXEXP.MAR; 1 
$f ; verified by M. Payne using her program ROOT16 on PDP-10. 
4 125 :TABDB: .WORD *0040200,0 ; 2*#(0/16) = 1.0 
4 150 ; «WORD OQ 
40 151; “WORD *604 305. 20125303 ; 2**(1/16) 
4 152 ; «WORD *006 +200 4) 
4 155 ; «WORD *004 § oO he! 1 3 2**(2/16) 
4 154 ; «WORD *0161752,°0105727 
4 155 ; «WORD “004 “014617 ; 3 2e0(3/16) 
4 1 § ; «WORD “007 eet 
4 157 ; -WORD %*004 "0 289 3 2%*(4/16) 
4 1 8 3 «WORD “005 ,“01 3¢ 1 
4 159 ; «WORD *004 9.20 § 3 20*(5/16) 
4 160 ; «WORD *006 oa0le & 
4 161 ; «WORD %*0040245,°0177326 3 2%*(6/16) 
04 162 ; «WORD *0124 81 + 00308 1 
04 163 ; «WORD antes 5,*005407 3: 2e*(7/16) 
04 164 ; «WORD “*016 oe: 1 
ot 165 ; «WORD 2002088 "000236 3; 2**(8/16) 
04 196 ; «WORD %*0031771,°015714 
040 167; «WORD apyes ENS ot ot 3 2%*(9/16) 
0040 168 ; «WORD “O11 Be + 0806 7 
0040 192 3 «WORD *0040305,*°00634 § : 2*2(10/16) 
0040 170; «WORD *0010 $2.08 33 
Once 171 ; -WORD *004 8.00 14 3; 20*(11/16) 
04 17 3 «WORD *0012437,°O1 9 
0040 173; «WORD %0040327,°0042374 : 2%*(12/16) 
0040 174; «WORD *0145526,°0116553 
0040 175; «WORD %0040340,°0146336 : 2%0(13/16) 
0040 176 3 «WORD *0166052,°011 $04 
Bpep 177 ; «WORD apres 20140306 3; 2%*(14/16) 
04 178 ; «WORD “016 ,“00 gf) 
Beh 179 ; «WORD *0040365,*0022575 3; 20*(15/16) 
4 180 ; «WORD *0012444,°0103314 
O40 183 
et 188 ; Table to be used for scaling. These constants are truncated 
40 184 ; conversions of the D_floating constants above to G_floating 
040 185 ; values. 
ae 
00000000 0000% 19 04 1 3 TABGB: .QUAD %*x0000000 00004018 : 2%*(0/16) = 1.0 
9708 CF9 401 4 1 -QUAD “XxX gor gee 401 3 2*2(1/16) 
17A 74 g5ke 11 : 190 -QUAD “*X517A 7p $557 11 3 wet TAY fs 
tr ¥4 A401 191 + QUAD alse ace A401 : 2e*(3/16) 
150A31 06FE401 6 136 QUAD °*XB7150A31 Ant} > 2**(4/16) 
4 +9 DEA6401 : 19 QUAD “X34 na DEA6401 3; 2**(5/16) 
A27D BFDA4O14 194 -QUAD “X2A27D FDASO14 3 2**(6/16) 
429004 ryt Eo 15 7 195 -QUAD *X5429004 74015 3 2*0(7/16) 
cf 7F AOSES 1 138 -QUAD *X aces TFA fe 6 ; eetaeizs 
1 B A114401 19 -QUAD “Xx BA114401 > 2%*(9/16) 
AODB422A are? 3 198 -QUAD “XA 44 Race ¢ 13 > 2**(10/16) 
FO90B2A3 (491401 9 19 QUAD “XF 3 A3C491401 > 2*0(11/16) 
3 AD99SA 589r 1A A ? QUAD “XDS35AD ators 1A 3 ett att 
9CDDB5 401C A QUAD “Xx 2500 5 401C > 2**(15/16) 
A487DCFB 58184010 § QUAD “XA4870CFB5818401D > 2**(14/16) 
QOD9AZA4 AGAF4OIE QUAD “X9OD9AZAGAGAF4OIE : 2e*(15/16) 
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& 16-SEP-1984 04:59 
T ; Declarative Part of Modu g-3 ety y 94:95:82 
Ms 
: $ 3 Polynomial coefficient tables for POLYF 
C : EXPTAB 
27¢e 3, C «WORD *0 243 "01146444 s 
F ¢ 19 «WORD *0O Ye 030081 838 ; 
FDF c8 1 » WORD *0 ; 65,°01 144 : 
721 1 35 \§ «WORD %*0037061,°00710 3 
0 0 D 1 «WORD 5 
000000 D4 18 EXPLEN = <.-EXPTAB>/4 
D4 18 EXPTB1 
HEH in 0D4 1 «WORD *0037052,°0131463 $ 
B5 F2A 00D8 18 «WORD “0 Iie ,“01 $3 5 F 
FFFFE SFFF 45 1 «WORD *0037777,°0177777 : 
FFFF 4O7F OOE 0 WORD *0040177,°0177777 : 
0000 anee OES 1 -WORD *0040200,0 g 
0000000 ors § EXPLN1 = <.-EXPTB1>/4 
OE 4 F_16L0G2_E: $ 
AA3B 4288 O0E8 5 ~WORD %*0041270,°0125073 
OEC $ F_LN2_OV_16_HI: ; 
70003E31 OOEC ~LONG %*x70003E31 
ta 8 F_LN2_0V_16_LO: ; 
FDF43705 OOF 9 ~ LONG “XFDF43705 
OOFS4 0 
OOF 4 1 
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1.00000000254251 


; LOG2(E) * 16 
; High 13 bits |n2/16 
; Low bits of ln2/16 
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-SBTTL MTHSEXP = Standard Single Precision Floating EXP 


Ww 


++ 
FUNCTIONAL DESCRIPTION: 
; EXP = single precision floating point function 


7 bit accuracy. 


: CALLING SEQUENCE: 
Exponential.wf.v = MTHSEXP(x.rf.r) 
INPUT PARAMETERS: 


LONG = 4 ; define songeere multiplier 
x = 1 * LONG ; Contents of x is the argument 


; IMPLICIT INPUTS: none 
; OUTPUT PARAMETERS: 
VALUE: floating exponential of the argument 
IMPLICIT OUTPUTS: none 
SIDE EFFECTS: 
Signals: MTHS_FLOOVEMAT if X > 88.028 with reserved operand in RO/R1 (copied 


: haath is computed using Chebyshev approximation 1001: about 


00000004 
00000004 
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65 ; to the signal mechanism vector CHFSL_MCH_RO/R1 by LIBSSIGNAL). Associated 
6 message is: ‘FLOATING OVERFLOW IN MATH LIBRARY’. Result is reserved operand 
6 -0.0 unless a user puge! red (or any) error handler changes CHFSL_MCH_RO/R1. 
68 MTHS_FLOUNDMAT if X =< -89.416 and caller has hardware enable Set. 
9 ; The result is set to +0.0. Associated message is: ‘FLOATING UNDERFLOW 
oy IN MATH LIBRARY’ 
7 NOTE: This procedure disables floating point underflow, enable integer 
7 overflow, causes no floating overflow or other arithmetic traps, an 
ee preserves enables across the call. 
7 == 
4 
401C 7 -ENTRY MTHSEXP, ACMASK ; standard call-by-reference entry 
80 ; disable DV (and FU), enable IV 
81 MTHSFLAG_JACKET ; flag that this is a jacket procedure 


oOooo 


6D  O00000000'GF 9E MOVAB G*MTHSSJACKET_HND, (FP) 
; set handler address to jacket 


handler 


Soooo 


; in case of an error in special JSB 
; routine 
; RO = user's arg 
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MOVF @x(AP), RO 
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Ex 
rec 


BSBB 
RET 


, 


G 19 
onential 1 
sion Floa 


MTHSEXP_RG 


6-88 


P=] 
P=} 


Be 99:88:32 


; RO = special 


AX/VMS Macro v04-00 
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UVXEXP.MAR; 1 
XP(RO) 


3; return - result in RO 
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; Sin 


ATHSERP. 


-v 
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4 51 4A 
51. 54S 
52. 51 COAF 4 
0 52 42 
51 BC AF 44 
0 51 4 
50 AD AF 4 
07. «11 


50 54 S50 29 AT AF 54 


FF75 CF 04 50 55 
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«SBTTL MTHSEXP_R4 = Special EXP 


routine 


; Special EXP - used by the standard, and direct interfaces. 


MTHSEXP_R4:: 
OVUND: bieus #*x8000, RO, R 
pee 


CALLING SEQUENCE: 
pove anything needed in RO:R4 


JSB MTHSEXP_RG 
return with result in RO 


input in RO 


Note: This routine is written to avoid causing any integer overflows, 
floating overflows, or floating underflows or divide by 0 conditions, 


whether enabled or not. 
REGISTERS USED: 


RO = floating argument, then result 


R2 = diddled exponent 
R5 = scratch 
R4 = integer part of X * LOG2(E)* 


#*X3E00, R2, R 
#*xX5B0, R3 
BLSSU = SMTST 


CMPW R2, #*X4280 ; 
BLSS EVAL $ 


16 


special EXP routine 

Preliminary test for over/underflow 
R3 = (4+exponent) + Ist 7 fract bits 
Compare iX: with 88 

to more tests if LSSU 

else, -4 < unbiased exp < 8 

no exceptions in EMODF or POLYF 


Check for loss of significance in 
EMOD ( iX! >= 284 ) 
No loss of significance 


iXi >= 2e#4, EMOD will lose significance so the interger and fractional 
parts of X*16/ln2 must be obtained in seperate steps. 


MULF3 F_16L0G2_E, RO, R1 
CVTFL «RT, RG 


CVTILF R4, RI : 
MULES F_LN2_OV_16_HI, R1, R2 ; 
MULF = F_LN2_OV_16_LO, R1 : 
SUBF RT, R : 
MULF F_{6L0G2_E, RO ; 
BRB POLY : 


EVAL: EMODF F_16L0G2_E, #X_51, RO, R4 


POLY: POLYF RO, #EXPLEN-1, EXPTAB 


BICL3 #-16, R4, R2 
MULF - TABLECR24, RO 


Get yopenee part of X*16/Ln2 in 
R4 (=]+J) as a longword and in 
R1 in F format 

Get fraction part of X*16/Ln2 = 
J6/inest KX = (14J)*ln2/16 J 
n . 


, RO 
get X*16*LG2(E) with 
pateper part in R4 (=]+4J) 
fraction = W in RO/R1 
evaluate polynomial ap- 
roximation with POLY. 
coefficients. 


R2 = J 
else MUL by 2**(J/16) 
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MTHSEXP_RG = Special EXP routine 6-SEP=1 MTHRTL.SRCJUVXEXP.MAR; 1 (5) 
159 45 ; Replace the following two Lines of code with the corresponding 
! : r ; G_floating code. 
1 48 : DDD TABDBER2J, RO ; add in DP 2*(y/16) 
} 4 49 ; CVTDF RO, R 
50 50 99FD 0159 $1 ‘ CVTFG RO, R 
50 ob cr4 é ED 130 36 ADDG¢ TAGB R2], RO : add in DP 2*#(J/16) 
4 OF cA 168 354 BICL #15, R64 R42 : 
3 0168 55 BEQL RETURN : if 120, then done 
50 6044 7E 190 36 MOVAQ (RO)CR4I], RO : shift I to exp position and 
171 3 MUL by 2**] by exponent addition 
OO7F B&F 30 81 171 28 CMPW RO, #*X7F ; test for over/under flow 
1§ 017 5 BLEQ EXCEPT : see what exception is if neg or = 0 
05 \e 9 RETURN: RSB 3; otherwise return result in RO 
179 6¢ SMTST: 
1 19 179 6 BLSS TOOBIG i exception if exp+4 > 11 
3400 8F 5 81 HE: 64 CMPW R2, #*x3400 ; eliminate underflow if exp <-24 
19 018 65 BLSS ; bypass if E**ARG = 1 
FFGC CF 04 50 55 018 06 POLYF RO, #EXPLN1I-1, EXPTB1 ; evaluate alternate polynomial 
05 Hy rt RSB 
50 08 50 0189 69 ONE: MOVF $*#1.0, RO ; E**ARG =1; store it 
05 018C 70 RSB 3 and return 
018D 71 
18> 37 
180 if : Handlers for software detected over/underflow conditions follow 
50 53 b18p 6 T00B1G: TSTF RO : if big ARG > 0 goto OVERFLOW 
28 18 Hb 4 BGEQ OVER 
a 23 : Underflow; if user has FU set, signal error. Always return 0.0 
0191 1 UNDER: 
52. DC «(0191 Be MOVPSL R2 3 Re = user's or jacket routine's PSL 
00000000'GF 00 FB 0193 8 CALLS #0, G*MTHSSJACKET_TST ; RO = TRUE if JSB from jacket routine 
04 50 9 OISA 4 BLBC R 3; branch if user did JS 
52 04 AD C 19D 5 MOVZWL SF$w_SAVE_PSW(FP), R2 : ast user PSL saved by CALL 
50 D4 O1AI $ 10$: CLRL RO ; RO = result. LIBSSIGNAL will save in 
1A ; CHFSL_MCH_RO/R1 so any handler can fixup 
0D 52 06 £1 OTA 8 BBC #6, R2, 208 ; has user enabled floating underflow? 
DD O1A PUSHL p > yes, return PC from special routine 
7E OO°8F QA OAD 90 MOVZBL #MTHSK_FLOUNDMAT, -(SP) ; trap code for hardware floa ing under f low 
1AD 91 3; convert to MTH$_FLOUNDMAT (32-bit VAX-11 
1AD 35 3 exception code) 
00000000'GF 02 Fe 1AD 9 CALLS #2, G*MTHSSSIGNAL z signal (condition, PC) 
0 18 ae 20$: RSB 3; return 
1B 38 EXCEPT: 
54 D5 O18 9 TSTL RG ; test sign of I; if 1 <0 
d&8 19 44 : BLSS UNDER ; go to underflow handler 
H+ 4 780 ; Signal floating overflow, return reserved operand, -0.0 


a “aay 


le Precision Floating Ex ouaties 16-SEP-1984 4:59 VAX/VMS Macro V04-00 Pp 1 Vv) 
P_R4 = Special EXP rout t he pease 7 tf 98 754 CMTHRTL.SRCJUVXEXP.MAR; 1 — (3) =( 
89 else process for overflow 

89 return PC from special routine 

8B uy | nO ores feeet ine over flow 


= — ; = reserved eperenc 
0. RO will be copied t 
signal mechanism vector CCHF SL _MCH_RO/R1) 
so can be fixed up by any error 
handler 
signal (condition, PC) 
return - RO restored from CHFSL_MCH_RO/R1 


00000000°GF 02 fe tS #2, G*MTHSSSIGNAL 


1 

E dD Q1 

a ae 
50. (0 ; : 
1 

1 

1 

1 

1 

1 

8 

1 


i 
$ 
OVER 

PUSHL (SP) 
MOVZBL #MTHSK_FLOOVEMAT, -(SP) 
ASHQ aT, R 

1 

1 

1 

1 

1 


4 
6 
404 
4 
4 
4 
4 
4 
4 
4 
6 
4 
4 


AOOOOOOOOW 


DOW > Wiwiw m 


FUN - OOONAUSE 


«END 


UVXSEXP 
Symbol table 


Symbol table output 

Psect synopsis output 
Cross-reference output 
Assembler run totals 


K 1 
; Single Precision Floating E ccoatial 16-SEP-1984 :04:59 VAX/VMS M v04 “y 
’ wiles pa 7 960482 MTHRTL.S RE SovX VXEXP.MAR; 1 


ACMASK = 
AL R os 
EXCEPT R 
EXPLEN = 
EXPLNi = 
EXPTAB R 
EXPTB1 g 
16 062 & 
F~LN2_OV_16_H1 # 
F~LN2_OV_16_LO a 0 
= 
MTIHSSJACKET_HND teeeenee 02 
MTHSSJACKET_TST teeeeenre 
MTHSSS1 L reerenee x 
THSEX 44 RG 
MTHSEXP_RG 000104 RG 
MTHSK_FCOOVEMAT eeeeeeee =X 
MTHS$K_FLOUND aeereeee =X 0 
ONE 0000189 R 0 
OVER 000001B9 R 0 
OVUND 00000104 R 0 
POLY 0000145 R 0 
RETURN 000178 R 0 
SFSW_SAVE_PSW = sig 
SMTST 00000179 R 0 
TABGB 00000040 R 0 
TABLE 00000000 R 0 
TOOBIG 0000018D R 0 
UNDER 00000191 R 02 
= 00000004 
x_51 = 00000029 
pom eemwmoe wea secre + 
: :_Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
ABS 00000000 < 8 00 -) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
S$ 00000000 < 0.) 4 -) NOPIC USR CON ABS LCL NOSHR EXE RD WRT NOVEC BYTE 
MTHSCODE 000001¢CB ¢« 45 2 . PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
eeusdusesseesesesooesses} 
H ; Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 3 0.10 0:0 “ 
zomnend processing 1; : . 0 “ 
Symbol table sort 
ass 2 
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uynte 
VAX-11 Macro Run Statistics 


The working set Limit was 1050 pages. 
reve bytes ‘ 6 pages) of virtual memory were used 7, buffer 
Ibere were pages of symbol gente, space ducing 18 of hold 
3 source l nes were read in Pass prod ob 
9 pages of virtual memory were ot, to define @ macros. 


Macro Library name Macros defined 
_$255$DUA28: CSYSLIBJSTARLET.MLB; 2 4 
88 GETS were required to define 4 macros. 


There were no errors, warnings or information messages. 


L 10 
; Single Precision Floating Exponential 04-00 
RCI UVXEXP.MAR; 1 


he intermediate 
9 non-local and 


gre 
local symbols. 
ect records in Pass 


cme me meee een aenan ane m ncaa es + 


! Macro Library statistics ; 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:UVXEXP/OBJ=OBJ$:UVXEXP MSRC$:MTHJACKET/UPDATE=(ENH$ :MTHJACKET) #MSRC$:UV 


| 
———eee 
LL 


| 0265 AH-BT13A-SE ) NT CORPORATION 
gave VAX/VMS V4.0 D PROPRIETARY 


